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Introduction 


In this paper, we will present and apply a methodology for identi- 
fying from experimental data the parameters in a multi -input/multi -output 
model of manual control, suitable for analysis of complex tasks. In this 
context, "complex tasks" refers to tasks in which multiple loop closures 
are expected to be present, such as multi -axis tracking or aircraft 
landing approach, as opposed to compensatory tracking in the laboratory, 
for example. The structure of the model is compatible with the well- 
known^^ optimal -control model (OCM) of the human operator, and among the 
model -related parameters we seek to obtain is an estimate of the manual 
controller's objective function "weightings". This is considered important 
because by doing so, the modeler may obtain insight into the operator's 
strategy and perception of the task. Furthermore, the magnitude this 
function takes on has been hypothesized 1 - * J to correlate with the 
operator's subjective rating of this task. Hence, an experimentally 
determined metric related to the subjective assessment of the task may 
hopeful ly result. 

In addition to establishing a model structure useful for identifi- 
cation, we will also evaluate two procedures for the determination of 
the desired model parameters. One uses frequency-domain measurements, 
and as a result is similar to previous methods. However, a variation 
on this technique will be presented to facilitate multi -input/multi -output 
model determination. The second procedure- proposed is based entirely 
on time-domain data. As a result, the constraints on the experimental 
procedure, such as special tracking signals, required in the frequency- 
domain approaches are avoided. Results from both methods will be presented 
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Model Structure 


Since the model structure to be used is to be compatible with the 
OCM, we will briefly note its key features. Readers unfamiliar with this 
modeling approach are referred to the references. The hypothesis upon 
which the model is based is that the well trained , well motivated human 
controller chooses his control inputs (e.g. stick force) to meet his 
(internal) objective in the task, subject to his human limitations. This 
objective is further assumed to be expressible in terms of a quadratic 


"cost" function 
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Q,F,R 


= vector of human's observed variables (e.g 
eration) 

= vector of his control inputs 
= Controller-Selected (internal) weightings 


attitude, accel- 


The human limitations modeled include information-acquisition and 
processing time delay, observation and control input errors, and neuro- 
muscular dynamics. A block diagram of the resulting model structure is 
shown in Figure 1. 

The components of this model may be grouped into two parts, one 
dealing with the information acquisition' and state estimation, and one 
related to the control law or control policy operating on the estimated 
state. As has been shown in the references on this modeling approach, 
the "solution" for the human's control inputs, as predicted by the 
model, is expressed as 
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Figure 1. Optimal Control Model 




( 1 ) 


u p = - G x x ■ G u u p + v u 

A 

where x = internal estimate of the system states 
G ,G = control gain matrices 

X u 

v u = motor noise, or control input errors 
The system dynamics are taken as 


x = Ax + BUp ; y = Cx 
where y = vector of system response 

and if a tracking task is considered, the dynamics of the tracking signal 
vector y may be represented as 

V 


y c V c + Dw 

where w is a disturbance input of "white" noise. Usually, the tracking 
dynamics are combined with the system dynamics, resulting in an augmented 
state vector col[x,y c ]. Finally, the manual controller is considered 
to observe delayed system responses and commands, with some observation 
error, or 


V 


C x(t-i) 
y c (t-T) . 


+ V. 


where v y = vector of observation errors. Usually, tracking error 
e = y - y is observed, plus the commands v themselves if the task is that 

V 

of pursuit.- Additionally, other system responses may be observed but 
not actually regulated or tracked. Therefore, for our purposes we will 
arrange the human's observation vector as follows (dropping the t-x 
here for brevity) 

yj - [J el yj, ,J] ♦ »J 
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with y Q = C Q x representing observed responses other than errors and 
commands. Clearly, the above expression can always be represented in the 
form 

y p = Cx(t - t) + v y 

T T 

where x = [y c , x] and C is partitioned accordingly to yield 



Reference 1, for example gives closed-form expressions for the 
state covariance matrix E{xx^} for this structure, under assumptions of 
independence and "whiteness" on w, v , and v . A compatible frequency 
domain representation of the manual controller may also be obtained that 
effectively has the following form 

u p (s) = r'ujHUKYpis) + y + t'^s)^ O) 

where Y (s) = Laplace Transform of y D (t) (not delayed) 

r r 

T~*(s) = Neuromotor Dynamics (= [G^s + I]~* if Eqn. 1 is considered) 
H(s) = Manual Controller Transfer Function Matrix (Refer to Fig. 1) 

N ,N = Noise Vectors - Related to v and v 
y ^ y u 

Also, the command and system dynamics expressed as 

Y c (s) = [si - A c ]~ 1 DW(s ) = * c (s)DW(s) 
and 

X ( s ) = [si - A] -1 BU p (s) = $(s)BU p (s ) 
may be combined to form 
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X(S) = 


\(sf 


** c (s)DW(s)' 

X(s) 


$(s)BU p (s) 


Then 


Y p (s) = Cx(s) 



( 4 ) 


where recall C is partitioned as in Equation 2. If we now partition 
H(s) and N to be compatible with Y (s), we may let 

y p 

H(s) = [H e (s), H c (s), H Q (s )] 


With this structure, we may represent the system as in Fig. 2, where we 
have used the following definitions 

C c => [I 4 .0] (I = Identity) 

c e => [i:-c r ] 

and . 

c o C° :C o ] 

to make the matrices have compatible dimensions. Note that not only is 
this structure consistent with the OCM with multiple inputs and outputs, 
but in the scalar case with H = H = 0 the structure also reduces to 

V w 

the conventional compensatory tracking block diagram. Shown in Fig. 3 
is this simpler case. 

For experimentally estimating a scalar Y (jio) ; measurements are taken 
of e c (t), e(t) , and u (t). Spectral analysis is then performed to 
obtain 
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Figure 2. General Model Structure 













$ (jw) = Cross spectrum between 0„ and u 
cu r c p 

4> ce (jco) = Cross spectrum between 0 C and e 
and the desired relation is 

Vp(jaj) = $ cu (ja))/$ Ce (jai) 

which can be derived from block diagrm algebra (see for example Refs. 

4 and 5). 

Now, as discussed in Ref. 5, special experimental conditions must 
be invoked to identify multiple human operator transfer functions, as 
in H(s) discussed above. Specifically, independent excitation of all 
inputs to H ( s ) must be present, and this is frequently not possible in 
many practical situations. However, some alternate expressions will be 
developed which yield identifiable transfer functions directly related 
to the general model structure discussed here, but are not human 
operator transfer functions, like Y p (jw) in the scalar case cited above. 
Referring to Eqns. 3 and 4 above, or equivalently Fig. 2, we have 

U (s) = T“ 1 (s)[H e (s)(e(s) + N £ ) + H c (s)(Y c (s) + N c ) 

+ H 0 (s)(Y q (s) + N 0 ) + N u ] 
e(s) = Y c (s) - C r *(s)BU p (s) 

Y q (s) = C 0 $(s)BU p (s) 

Substituting e(s) and Y (s) into the first expression, and solving for 
U (s) yields 
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U p (s) - CT uc (s)]Y c (s) + ,[T uy (s)]N y+ CT uu (s)]N u 

( 5 ) 

where = [n[, nJ, N q ] 

T uc (s) = [I + r 1 (s){H e (s)C r - H 0 (s)C o }$(s)B]“ 1 [r 1 (s){H £ (s)+H c (s)}] 

T uy (s) = [I + T~ 1 (s){H £ (s)C r - H 0 (s)C 0 }$(s)B]' 1 [T‘ 1 (s)[H e (s):H c (s):H 0 (s)]] 

T uu (s) = [I + T n 1(s){H e (s)C r ' H 0 ( s ) c 0 H(s)B] " lT n 1(s) 

Substitution back into the relation for e(s) and Y Q (s) yields 

e(s) - [T cC (s)]Y c (s) ♦ [T Ey (=)]N y + [T eU (s)3N u (6) 

Y 0 (s) = [T 0C (S)]Y C (S) + [T oy (s)]N y ♦ [T 0U (S)]N U 

where 

T £C (S) - [I - C r *(s)BT uc (s)] 

V $) - -C r »(s)BT uy (s) 

T eu (s) = -C r »(s)BT uu (s) 

T oc (s) - C 0 <f(s)BT uc (s) 
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and 


V s) = V< s > BT uo< s > 

T ou (s) = C o *(s)BT uu (s) 


Now the three transfer function matrices T , T , and T are related 

uc EC oc 

to the matrices of cross-spectra between Up and y c , e and y c> and y Q and 
y , respectively, assuming the "noises" N and N are uncorrelated with 
y Q . Or using matrix notation 


T uc^> = ^ y U y 

J c p ^c^c 

(?) 

T oc (jO)) = O y y (jl0)][* y y (U)]' 1 
0 0 0 0 

So, if frequency-domain data were used to estimate the above spectra, 

the transfer functions in Eqn. 7 may be identified, but not necessarily 

the elements of H(s). However, these identifiable transfer functions, 

due to their direct relationship to the OCM, for example, may be used 

for model identification and/or validation in exactly the same manner 

estimates for H(s) may be used, so they are just as meaningful. 

Additionally, referring back to Equation 5, under the assumption 

that the noise vectors N and N consist of elements mutually uncorrelated, 

%y ^ 

and uncorrelated with y , a model -related ‘expression for the power of 

V» 

the remnant in each of the i'th components of Up is expressible as 


*rr (u) = £l T u.y. (jw) | 2 yy. (a)) 

3 n ■Vj 


+ H T u u «-)l 2 * u u ^ 
K u i u k u k u k 


( 8 ) 
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where $ (w) = Power spectrum of the j'th element in the noise N 

y j y j y 

$ (w) = Power spectrum of the K'th element in the noise N . 

u k u k 

So if $^ r (a>) is estimated experimentally, it is relatable to the model- 
based parameters on the right hand side of the above equation for further 
model comparison or validation. Similar expressions for all the above 
development are available in Ref. 6, for further reference. 

Parameter Search Technique 

Now that the model structure is obtained to allow direct comparisons 

between measured variables and their model-based counterparts, attention 

is now turned to obtaining the parameter set of interest. This parameter 

set, denoted p, consists of the "independent" variables of the model, 

such as objective function weights Q and R, time delay x, and noise 

covariance matrices C and C , for example. We will make direct 

n y n u 

application of the quasi -Newton search approach of Refs. 5, 7 and 8, 
with two variations fundamental to our purpose. The first is that in 
the above references, a scalar objective function weight on tracking 
error alone was used exclusively, while we desire to estimate more complex 
expressions for the cost, or weighting matrices. Secondly, we will 
compare using two forms of experimental data, one strictly time domain 
and the other frequency domain, to determine if using only time domain 
data leads to sufficiently accurate results. This is desirable since a 
purely time domain approach is simpler and greatly reduces the requirement 
on the experimental technique for obtaining the required data. 

The scheme is implemented to minimize a scalar matching cost of 
the form 
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X Vi 

where e. is the difference between the i n measured data point and the 
corresponding model prediction, w. is a weighting coefficient. Or in 
matrix form: 

M = e T We 

with e = colCe^, ...], W = diagfWj], 

For a trial set of model parameters p^, we have its corresponding 
modeling cost 

Mj = ejkle^ 

For a new set of parameters Pg = p^ + Ap, we obtain a new modeling 
error Ae, related to Ap by 

Ae = QAp 
9e, 

where q(i, j) = — - can be obtained by a numerical perturbation of the 

3 Pj 

model. The change in the parameter vector Ap yielding the minimum 
modeling error, given the initial vector e^ and the assumption of 
linearity between aM and Ap is 

Ap = -[Q T WQ]“ 1 Q T We 1 

Thus an iteration procedure is established, which proceeds until no more 
improvement in matching cost M, or the required changes in the parameters 
in Ap are very small. 

In addition to obtaining the best match to a given set of data, we 
also wish to determine some measure of the reliability of the identified 


129 



parameter values. A qualitative indication of parameter estimation relia- 
bility can often be obtained through sensitivity analysis relating changes 
in the scalar matching cost to perturbations in the model parameters. In 
general, estimates of parameters that have a high impact on the matching 
cost can be considered more reliable than estimates of parameters having 
a smaller impact. 

As shown in Ref. 8, this sensitivity may be estimated from the relation 
aM-j = V T Q T WQV(Ap.) 2 

X.L 

where V is a column vector that has a value of unity for the i n element 
and values for remaining elements V r as determined from 

V r - -CQjwQ r 3 _1 qj[wq 1 

where q.. = col[q i j, q i 2 , ...] and the subscript r indicates vectors and 
matrices which omission of the ith row and column. 

Pursuit Tracking Analysis 

For comparing the time-vs. frequency -domain data for model determin- 
ation, and to relate the above methodology to an established situation, 
a single axis pursuit tracking task is considered. ^ Subjects tracked a 
command signal generated by a sum of sinusoids 

10 

© r = l A. si n(u) . t + <f>.) . 
c i=1 i i 1 

for 100 seconds, with the frequencies w. evenly spaced between 0.25-17 
rad/sec, and amplitudes A^ selected such that the spectrum of the 
command approximated a random signal generated by 

0 /w = 1 

c s^ + 3s + 2.25 
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2 

with "white" noise w intensity taken to be a = 13.5 6(t), - s(t) 
a delta function. 

In addition to 0 , the subjects observed the plant response e(t), 

and therefore the error e - 0, where the two olants ( 0/6 ) were K/s 

2 P 
and K/s . A representative block diagram is shown in Figure 4. (Note 

the correspondence between this block' diagram and that of Figure 2.) 

Since error, e , and 0 are not all linearly independent, only two need 

be included for observation. Therefore, the subjects observation vector 

may be taken as 

yj - [s. c, e c , e c ] 

for both K/s and K/s plants. Finally, including the subjects' control 
input 6 in the state vector x, we may define 

r 

For K/s , x T = [0 C , 0 C , 0 , 5 ] 

For K/s 2 , x T = [0 C , 0 C , 0 , 0, 6 ] 

Referring back to Eqns. 5 and 6, one may consider T (s) and T (s) to 

uc ec 

be scalars, 

T uc (s) = [H 1 + sH 2 + H 3 + sH 4 ][1 + F (Hi+sH 2 )/(T n s + 1 )]" 1 

n 

T ec (s) = [1 - F(H 3 +sH 4 )/(T n s+l)][Jl + F(H 1 +sH 2 )/(T n s+l)]“ 1 

From the experimental data, the state covariance matrices E[xx T ] 

were estimated, as well as the cross-spectra between 0„ and 6. and 0 

c p c 

and e (error), or g (jw) and e (jw). Finally, although not 

w c°p w c e 

possible in more complex situations, since T (s) and T (s) are scalars 

uc ec 

in this case, an effective operator transfer function may be defined as 
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y P eff (dw) = T uc (joj)/T e c (ja)) 


= $ 0 & (ja))/^ e (jw) 
c p 


V 


These time and frequency-domain results were used for the quasi-Newton 
parameter search to estimate 

P T ■ [q E . q-. x n , t, Cj . C ] 

J ” 

where q , q- = objective function weights on error and error rate 
= neuromotor time constant 
t = observation time delay 

C p = observation noise intensities (expressed as noise-to- 
signal ratios in dB relative to the variance of each 
observation) 

C = motor noise intensity (expressed as noise-to-signal ratio 
u 

in dB relative to control input variance) 

Two separate parameter searches were performed. One used only the 
state covariance matrix for computing model matching cost, or 


M t = 


l N 2 

w 

1J 


where X.. 

* J 

A 

X ij 


a . . 


element in experimentally-obtained state covariance matrix 
corresponding element from the covariance matrix from the 
model 

standard deviation in the experimental value of X.. over 

J 

the repeated runs. 
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The second used only the freauency-domain result for y (ju>. ), 

p eff 1 


N 


G.-G. 2 2 R.-R. 2 

r 1 j ./1 l \ j /1 


M r = ¥ I K-*— 1 ) + (- 1 — L ) + (— — -) ] 
F N 1 "G, °r 


where G. ip. = |y (jay. )|, and arg y_ (ju>. ),• <*>• innut freouencies in 
u 1 p eff 1 p eff 1 1 

command signal, measured experimentally from spectra 

A A 

<p= corresponding magnitude and phase of the model -estimated 
transfer function 

R.j = estimated power of the remnant in the control input fip, 

from experiment. Obtained from the spectrum of S at 

P 

frequencies other than those in the command. 

A 

R.j = remnant power obtained from the model (or Eqn. 8) 
a G. ,C V ,a R. = standard deviation of the experimental data 

The estimates for desired model parameters p obtained using both 
approaches are listed in Tables 1 and 2, for the K/s and K/s 2 plants, 
respectively. Note the estimated values of the parameters do not differ 
significantly between the results obtained from minimizing M T (time domain) 
and those from minimizing Mp (frequency domain). In some cases, the 
sensitivities in these costs to small relative changes in these parameters 
do vary, depending on whether frequency or time domain data is used. 

Another interesting result is the comparison between the state 

covariance matrices obtained from the frequency-data - matched model and 

the time-data - matched model. The results for the K/s plant are given 

2 

in Table 3, while those for K/s are shown in Table 4. These results show 
not only excellent agreement with simulation results, but the result 
from the frequency-domain match agrees very well with the time domain 
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Time Match 


Frequency Match 


QCM p i lot— re 1 a ted 
parameters 

identification 

results 

sensitivity 

identification 





Results 

time delay, T 

. 14 

sec 

. 47 

. 13 

sec 

weighting on error, q^ 

2707. 


. 19 

2380. 


weighting on error, q» 
rate e 

348. 


. 09 

377. 


motor noise , c 

U 

-18. 

db 

1. 39 

-18. 

db 

observation noise 
of command, c 

^c 

-6. 

db 

. 21 

-7. 

db 

observation noise 
of command rate* C* 

^c 

6 

t-4 

1 

db 

. 68 

-10. 

db 

observation noise 
of error , c 

e 

-11. 

db 

. 02 

-11. 

db 

observation noise 
of error rate, c. 

-10. 

db 

. 05 

-10. 

db 


neuromotor lag, x n .09 sec 


sensitivity 

10 . 

. 2 
4. 1 

3. 2 
2 . 2 

2. 9 

3. 5 
3. 4 


Table 1. Matching Results - K/s 


Plant 
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Time Match 


Frequency Match 


OCM pilot-related 
parameters 

identification s 

results 

ensitivity 

identification 

result 

sensitivity 

time delay, x 

.10 

sec 

. 65 

. 09 sec 


10. 

weighting on error, q 

E 

2130. 


. 25 

2320. 


. 1 

weighting on error 

319. 


. 18 

380. 


2. 8 

rate , q* 





motor noise, c 

U 

U 

-19. 

db 

1. 20 

-19. db 


. 3 

observation noise 
on command, C 

-7. 1 

db 

. 48 

-5.8 db 


2. 1 

observation noise 
on command rate, c* 

-7. 5 

db 

. 34 

-9.5 db 


1. 4 

observation noise 
on error , c 

e 

-15. 

db 

. 67 

-11. db 


1. 

observation noise 

on error rate , C* 

e 

—9. 7 

db 

. 55 

-10. db 


1. 5 

neuromotor lag 

.20 sec 






Table 2. 

Matching Results - K/s z Plant 






Table 3. Augmented State Covariance Matrix for K/s Plant 


simulation result: 
1. 0 

0. 

. 80 

. 46 

0 c (deg) 

O. 

2. 25 

46 

. 71 

© c (deg/sec) 

. 80 

46 

. 93 

1 

G (deg) 

. 46 

. 71 

1 

3. 7 

« p (in) 

frequency do/na in match : 
1.0 0. 

. 77 

. 28 


0. 

2. 25 

28 

. 91 


. 77 

28 

. 76 

0. 


• 28 

. 91 

. 0 

3. 5 


time domain match: 
1. 0 

0. 

. 79 

. 29 


0. 

2. 25 

29 

. 91 


, 79 

29 

. 80 

0. 


. 29 

. 91 

. 0 

4. 0 
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Table 4. Augmented State Covariance Matrix for K/s 2 Plant 


simulation result: 


1. 0 

O. 

. 73 

. 46 

3 

© c (deg) 

0. 

2. 25 

48 

. 13 

1. 8 

0 c (deg/sec) 

. 73 

48 

1. 27 

0. 

-3. 

0 (deg) 

. 46 

13 

0. 

2. 8 

32 

0 (deg/sec) 

3 

1. 8 

-3. 

-. 32 

35. 

V ©"> 

frequency domain 

match : 




1. 0 

0. 

. 67 

. 39 

-. 36 


0. 

2. 25 

39 

. 36 

2. O 


. 67 

39 

1. 18 

0. 

-2. 7 


. 39 

. 36 

0. 

2. 7 

0. 


-. 36 

2. O 

-2. 7 

0. 

38. 


time d 

omain match: 





1. 0 

o. 

. 73 

. 44 

-. 32 


0. 

2. 25 

44 

.‘32 

1. 9 


. 73 

-. 44 

1. 24 

0. 

-3. 0 


. 44 

. 32 

0. 

2. 9 

0. 


32 

1. 9 

-3. 0 

0. 

41. 
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model, obtained by matching these statistics. 

On the other hand, the frequency -matched model, as expected, matches 
that experimental data well, as shown in Figures 5 and 6. Note, further- 
more, that the time-matched model does not do a poor job of matching 
this data as well . 

From the above results, the following is noted: 

1. The model obtained from time-domain matching is very close to the 
model obtained using frequency-domain data. 

2. The sensitivity of the match to model parameter variations, however, 

differs between the time-and frequency-domain matches. 

3. From the frequency-domain matches especially, the sensitivity of 
the match to variation in the cost function weighting on error 
rate, q^, is quite large. This indicates that including this 
parameter in the cost function is significant. 


Multi -Axis Tracking Analysis 


As a final example, we will summarize the results of an analysis 

rgl 

r-omniav ni. t.ci/ *- J The task involves fixed-base 


of a complex multi -axis tracking task. 


simulated air-to-air tracking, with the display symbology as shown in 
Fig. 7. The sight symbol (box) is dynamic, representing a lead-computing 
sight. It's position relative to the fixed screen reference is defined 


by the coordinates x^ and x^. The relative position of the target 
is defined by and A nd the relative bank angle between 
the target and attacker is indicated by the target's bank angle on the 


screen. (Note, <j> re -j is zero for situation shown in the figure.) The 
linearized system dynamics are representative of tracking during a 4g, 


constant altitude turn. The input (or command) driving the closed-loop 
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CRT Display for Multiaxis Air-to-Air 
Tracking Task 

Fig. 7 ' 
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system is the target's inertial (not relative) bank angle <j>y, which is 
generated by the relation 

4>y = -1 /t (f>y + w 


with t = 13 sec., and the intensity of the random w selected to yield 
an rms value of <j>y of 5.25 degrees. 

One selected set of the pilot's observed variables is 




EL’ e EL’ e AZ ! 


AZ* P EL’ P EL’ P AZ 




^AZ’ ^ 


where - x^, tracking errors 

<f> = attacker's bank angle 

Other combinations of observations could also be selected, and this set may 
not be optimum. Variations on this are under investigation. The pilot's 
control input is the stick and rudder, or fi E , 6 A , and s R . 

The model parameters will be estimated by a time-domain matching 
of the (16 x 16) state covariance matrix, including the three control 
inputs, obtained from several simulation runs. The parameter set to be 
discussed includes'the (3 x 3) T^ matrix (or G" 1 ) associated with the 
three control inputs, the cost function weights 








] 


and the noise intensities 


(c r , c_ , c. ) = Variances on motor noises 

6 e 6 a 6 r 

(c , c* , c , c* ) = Variances on measurement noises 
e El e El £ AZ e Az 


The variances on the noises associated with the additional measurements 
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were fixed at - 13 dB after some initial studies. 

As with the selected observation vector, the selection of cost 

function weights is based on subjective judgement, and one set may in 

fact be more meaningful than the other. For example, the use of a 

weighting on relative bank angle between target and attacker, rather 

than on g j and g^ could be considered. For the set selected here, 

however, the results are given in Table 5, and the T matrix is 

o n 


T n 


27 0 0 

0 .31 -.15 

0 -.15 .30 


(sec) 


for 


Up = 


6 a , 


Sr). 


Note the relatively high sensitivity on the cost weightings on g £L 
and g A2 in Table 5. This is consistent with the results of Harvey^ 10 ^ 
in his evaluation of a similar single-axis task, in that weightings on 
observations in addition to tracking error and error rate were significant 
in obtaining a good model match. This fact is basic to the desire to be 
able to identify more complex cost functions, as noted in the introduction. 

Finally, although this match used the simple-to'-obtain state 
covariance matrix, comparisons or matching of frequency domain data is 
certainly possible if available from the experiment. If not, the frequency 
domain results from the model is available as a "prediction" of those human 
operator characteristics. 

Note that slightly more general expressions for the transfer function 
matrices T uc (s), T £C (s) ancl T oc ( s ) in E( l ua ' t1 ' ons 5 and 6 result in the 


above example 


[ 6 ] 


This arises due to the fact that the system dynamics 
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Table 5 


Identification Result - Time Domain 
for Multiaxis Air— to-Air Tracking Task 


□CM p i lot— related identification sensitivity 


parameters 

results 



time delay, x 

. 13 

sec 

. 2 

weighting on 
elevation error) q^j 

1501. 


3. 4 

weighting on 
eleva. error rate , 

340. 


, 1 

weighting on 
azimuth error, q^ 

1741. 


2. 4 

weighting on 
azimu. error rate, q^ 

320. 


. 1 

weight, on target 
elevation angle, q 

3 E1 

1573. 


. 2 

weight, on target 
eleva. angle rate , q; 

B E1 

248. 


1 . 

weight, on target 

1556. 


. 2 

azimuth angle, q 0 

b a 




weight, on target 

226. 


1 . 

azimu. angle rate , q' 

6 a 

elevator noise, c. 

E 

-21. 3 

db 

2. 0 

aileron noise, c„ 

6 A 

-20. 6 

db 

. 4 

rudder noise, c r 

R 

-19. 2 

db 

1. 3 

meas. noise on 

-12. 8 

db 

. 8 

eleva. error, c 




e E 




meas. noise on 
eleva. error rate, c* 

-13. 2 

db 

1. 4 

e E 




meas. noise on 
azimuth error, c 

-13. 3 

db 

1 . 1 

e A 




meas. noise on 

-13. 2 

db 

. 9 


azimu. error rate , c* 

e A 
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are not decoupled into command and plant dynamics, as assumed previously. 
As a result, the equation for Up(s) and e(s) (Eqns. 5 and 6) are developed 
from the relation 


x (s) = [si - A]" 1 [BU p (s) + DW ( s ) ] 

where 



A. A'~ 


o’ 


D 

A = 

c c 

II 

ICQ 


* D = 



0 A 


B 


0 


In the development of Eqns. 5 and 6, A' c in A was assumed zero. With this 
change, the development of the desired matrices proceeds directly, along 
with modifying Figure 2 accordingly. 


Summary and Conclusions 

An approach has been presented for identifying and/or validating 
multi -input/multi -output models for the manual controller in complex 
tracking tasks. In the more general case, the conventional human 
describing functions may not be directly identifiable, but measurable 
transfer matrices directly related to the. model were derived. In terms 
of model identification or validation, these transfer matrices are just 
as useful and meaningful as the conventional describing functions. 

Model -parameter identification using strictly time-domain data 
was demonstrated to yield excellent results for the single-axis pursuit 
task. The use of this approach avoids the necessity of obtaining 
frequency domain data, sometimes a practical constraint. However, shown 
in Ref. 11, time-series techniques may be used effectively to obtain 
frequency-domain representations directly compatible with the parameter 
identification method presented here. Furthermore, the time-series methods 
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would appear to circumvent several of the practical problems in obtaining 
frequency-domain representations - such as the necessity to be able to 
define the command signal characteristics. Therefore, model parameter 
estimation using frequency-domain representations are certainly of 
interest, and will remain useful. 

The results obtained from evaluation of a two-axis air-to-air tracking task 
with complex, high-order dynamics were briefly noted, primarily to 
demonstrate the type of analysis possible with this approach. 
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